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The quasiparticle (QP) band structures of both strainless and strained monolayer M0S2 are inves- 
tigated using more accurate many body perturbation GW theory and maximally localized Wannier 
functions (MLWFs) approach. For the strainless state, the partially self-consistent GWq (scGWo) 
calculation demonstrates that monolayer M0S2 is a direct semiconductor, consistent with recent ex- 
perimental observations. By solving the Bethe-Salpeter equation (BSE) including excitonic effects, 
the predicted optical gap magnitude is consistent with available experimental data. We further 
predict that the strong exciton binding energy is in the range of 0.50—0.40 eV and also is nearly 
unchanged with increasing strain. In addition, the scGWo calculations also predict that monolayer 
M0S2 maintains its direct characteristic under the tensile strain up to 1.41%. As the tensile strain 
increases further, the electron effective mass decreases, the covalent bonding weakens, and a direct 
to indirect gap transition occurs. The scGWo and BSE calculations are also performed on monolayer 
WS2, similar characteristics are predicted and WS2 possesses the lightest effective mass at the same 
strain among monolayers Mo(S,Se) and W(S,Se). The present calculation results suggest a viable 
route to tune the electronic properties of monolayer transition-metal dichalcogenides (TMDs) using 
strain engineering for potential applications in high performance electronic devices. 

PACS numbers: 73.22.-f, 71.20.Nr, 71.35.-y 



I. INTRODUCTION 



Bulk TMDs consisting of two dimensional (2D) sheets 
bound to each other through weak van der Waals forces 
have been studied extensively owing to their potential ap- 
plications in photocatalysis ffl and catalysis fi [|. M0S2, 
WS 2 , MoSc 2 , and WSc 2 are examples of such TMDs. 
Recently, their 2D monolayer counterparts were success- 
fully fabricated using micromechanical cleavage method 
j| . Since then, these monolayer materials have attracted 
significant attention p|-[l3|. 

For monolayer M0S2, a strong photoluminescencc (PL) 
peak at about 1.90 eV, together with peaks at about 1.90 
and 2.05 eV of the adsorption spectrum indicated that 
M0S2 undergoes an indirect to direct band gap transi- 
tion when its bulk or multilayers form is replaced by a 
monolayer Theoretical studies employed density 

functional theory (DFT) method also predicted mono- 
layer M0S2 to have a direct gap of 1.78 eV ||. It is 
known however that DFT does not describe excited state 
of solids reliably. Furthermore, an important character in 
low-dimensional systems is their strong exciton binding 
due to the weak screening compared to bulk cases. There- 
fore, the good band gap agreement between theoretical 
and experimental results for monolayer M0S2 may be a 
mere coincidence. DFT calculations also showed that the 
band gap of monolayer M0S2 was sensitively dependent 
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on strain [|10| |12| , and the direct to indirect band gap 
transition was predicted at tensile strain of 0.3% p0| , 
while others reported at 1% |13[| . Furthermore, shifts 
of PL peak for the monolayer M0S2 were also observed 
experimentally, which was attributed to the strain intro- 
duced by covered oxides |T^| . In order to address above 
questions, it is important and necessary to employ more 
accurate calculation method beyond DFT to investigate 
the electronic structures of strained monolayer M0S2. 

The most common method to circumvent drawback of 
DFT is the GW approximation [fL5| , in which, self-energy 
operator £ contains all the electron-electron exchange 
and correlation effects. The scGWq approach, in which 
only the orbitals and eigenvalues in G are iterated, while 
W is fixed to the initial DFT Wo, was shown to be more 
accurate in many cases to predict band gaps of solids 
[pi. The off-diagonal components of the self-energy £ 
should be included in scGWq calculations, since this in- 
clusion has been proved particularly useful for materials 
such as NiO and MnO jl7|. It is noted that £ within 
GW approximation is defined only on a uniform k mesh 
in Brillouin-zone, due to its non-locality. Therefore, un- 
like DFT band structure plot, the QP eigenvalues at ar- 
bitrary k points along high symmetry lines cannot be 
performed directly Started from the sc GWq calcu- 
lation, the QP band structure can be interpolated using 
MLWFs approach. This combination was demonstrated 
to be accurate and efficient for the sc GW band structure 
|]l8| . The G W results were shown to agree well with the 
photoemission data [ fl9| , while in order to reproduce the 
experimental adsorption spectra, the consideration of at- 
traction between quasi-electron and quasi-hole (on top 
of GW approximation) by solving BSE is indispensable 
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|L9[ , particularly for the low-dimensional systems with 
strong excitonic effect. The main goal of this study is 
to accurately predict the QP band structures and opti- 
cal spectra of monolayer M0S2 as a function of strain by 
adopting the DFT-scGWVBSE approach. 

Strain in monolayer M0S2 can be produced either by 
epitaxy on a substrate or by mechanical loading. It is 
well-known that strain can be used to tunc the electronic 
properties of materials. This is particularly important 
for two-dimensional materials, which can sustain a large 
tensile strain. In fact, shifts of PL peak observed ex- 
perimentally in monolayer M0S2 was attributed to strain 
| |L4| , and the magnetic properties of M0S2 nanoribbons 
could be tuned by applying strain . 

By adopting the aforementioned approach, we system- 
atically investigate how the electronic structures and op- 
tical properties of monolayer M0S2 evolve as a function of 
strain. Based on the more accurate band structures, the 
effective masses of carriers are calculated. In addition, 
this calculation approach is also employed to investigate 
other monolayer TMDs, that is, WS 2 , MoSe 2 , and WSc 2 . 



II. DETAILS OF CALCULATION 

Our DFT calculations were performed by adopting the 
generalized gradient approximation (GGA) of PBE func- 
tional [^lj for the exchange correlation potential and the 
projector augmented wave (PAW) |0 method as imple- 
mented in the Vienna ab initio simulation package [ p3[ . 
12 valence electrons are included for both Mo and W 
pscudopotcntials. The electron wave function was ex- 
panded in a plane wave basis set with an energy cutoff of 
600 eV. A vacuum slab more than 15 A is added in the di- 
rection normal to the nanosheet plane. For the Brillouin 
zone integration, a 12x12x1 T centered Monkhorst-Pack 
fc-point mesh is used. In the following GW QP calcula- 
tions, both single-shot Go Wo and more accurate partially 
scGWo calculations are performed. 184 empty conduc- 
tion bands are included. The energy cutoff for the re- 
sponse function is set to be 300 eV, the obtained band 
gap value is almost identical to the case of 400 eV. The 
convergence of our calculations has been checked care- 
fully. For the Wannier band structure interpolation, d 
orbitals of Mo (W) and p orbitals of S (Se) are chosen 
for initial projections. Our BSE spectrum calculations 
are carried out on top of scGWq. The six highest va- 
lence bands and the eight lowest conduction bands were 
included as basis for the excitonic state. BSE was solved 
using the Tamm-Dancoff approximation. Notice that the 
applied strain in the present study is all equibiaxial, un- 
less stated otherwise. 



III. RESULTS AND DISCUSSIONS 

We first analyze the density of states (DOS) for mono- 
layer M0S2. The d orbitals of Mo and p orbitals of S con- 



tribute most to the states around the band gap, similar 
to previous studies fici|-|l2f . Fig. |l| shows the projected d 
orbitals of Mo and p orbitals of S as well as the decom- 
posed d orbitals for monolayer M0S2 at the experimental 
lattice of 3.160 A for bulk MoS 2 fHJ (Fig. 1(a) and (b)) 
and under 3% tensile strain (Fig. 1(c) and (d)), respec- 
tively Based on the DOS, the d orbitals of Mo and p 
orbitals of S are chosen as the initial projections in the 
Wannier interpolated method. Fig. |] shows the identi- 
cal DFT band structures of monolayer M0S2 obtained by 
the non-selfconsistent (Non-Scf) calculation at fixed po- 
tential and Wannier interpolation method, respectively, 
confirming that our choice of the initial projections and 
inner window energy is appropriate. Based on the good 
results for monolayer M0S2, the same procedure is also 
employed for remaining monolayer TMDs. 



A. QP band structures of strained monolayer M0S2 

The QP band structures of monolayer M0S2 at four lat- 
tice constants of 3.160, 3.190 (the optimized value from 
the present work), 3.255, and 3.350 A are plotted in Fig. 
H, corresponding to 0%, 1%, 3%, and 6% tensile strains 
(with reference to 3.160 A), respectively. As shown in 
Fig. |^(a), the band structure obtained by DFT for strain- 
less M0S2 is a direct band gap semiconductor with a band 
gap energy of 1.78 eV, while the indirect band gap of 2.49 
eV is predicted by Go Wq , which is consistent well with 
Olsen's results |9|. Obviously this Go Wq indirect band 
gap is contrary to the PL observations The QP 

band structures predicted by our scGWq calculation show 
that M0S2 is a K to K direct band gap semiconductor 
with a band gap energy of 2.80 eV. This prediction is in 
excellent agreement with the recent calculation for M0S2 
at the experimental lattice using full-potential linearized 
muffin-tin-orbital method (FP-LMTO) Q, which pre- 
dicted & K to K direct band gap of 2.76 eV. 

It should be noted that in the 2D materials, the exci- 
tonic effect is strong due to the weak screening. Thus it is 
important to consider the attraction between the quasi- 
electron and quasi-hole by solving the BSE discussed be- 
low in order to make the predicted optical gap consistent 
with the optical spectra. 

Fig. |^(b) shows the band structure of monolayer M0S2 
at 3.190 A corresponding to 1% strain. The DFT re- 
sult predicts the monolayer M0S2 to be an indirect band 
gap with K to Y of 1.67 cV. Previous DFT studies also 
found that monolayer M0S2 already becomes an indirect 
semiconductor under a tensile strain of 1% |13 . After 



GW correction, both of the Go Wq and scGWq QP band 
structures show that M0S2 is still a direct semiconductor 
with K to K band gaps of 2.50 and 2.66 eV, respectively. 
As the strain increases, shown in Fig. ||(c) and ||(d), 
the DFT, G Wq and scGWq all predict monolayer MoS 2 
to be indirect. The calculated indirect band gaps from 
DFT, G Wo and scGWo arc 1.20 (0.63), 2.19 (1.56), and 
2.23 (1.59) for monolayer MoS 2 under strain of 3% (6%), 
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FIG. 1: Projected density of states of d orbitals of Mo (a and c) and p orbitals of S (b and d) and decomposed d orbital of Mo 
for monolayer M0S2 at lattice constants of 3.160 (a and b) and 3.255 A (c and d), respectively. The latter corresponds to 3% 
strain. 




(a)a=3.160A (b)a=3.190 



FIG. 2: DFT band structures of monolayer M0S2 at lattice 
constant of 3.160 A obtained by the Non-Scf calculation and 
Wannier interpolation method. 





(c) a=3.255 A (d) a=3.350 A 





FIG. 3: DFT, Go Wo, and scGWo QP band structures for 
monolayer M0S2 at lattice constants of (a) 3.160, (b) 3.190 
(the optimized lattice constant from this work), (c) 3.255, and 
(d) 3.350 A corresponding to 0%, 1%, 3%, and 6% tensile 
strain (with reference to 3.160 A), respectively. The Fermi 
level is set to be zero. 



respectively. As shown in Fig. |], the value of band gap 
decreases as the tensile strain increases, accompanying a 
shift of valence band maximum ( VBM) from K to T point 
and resulting in a direct to indirect band gap transition. 

To further understand the direct to indirect gap transi- 
tion and determine this transition point, the band gaps as 
a function of tensile strain are plotted in Fig. |[ Clearly 
our DFT and scGWo results have the same trends, and 
accord well with reported DFT |l3| (cyan triangle) and 
scGW |24| (green solid square) results, respectively. The 
inset in Fig. ^ is the energy difference between the va- 
lence band top (VBT) at K and T points, which deter- 
mines the direct to indirect gap transition for monolayer 
M0S2 under tensile strain. The dependence of the en- 
ergy difference on the tensile strain is almost linear both 
for DFT and scGWq results, and the estimated direct to 



indirect transition occurs at about 0.94% and 1.41% for 
DFT and scGWo, respectively, marked by vertical black 
and red lines in Fig. [| The former is consistent well 
with reported DFT prediction of 1% JuJ, while a smaller 
value of 0.3% was obtained in Rcf.|10j. The difference 
originates from the relative position of VBT at K and T 
points mainly due to the pseudopotcntial and all-electron 
methods adopted at DFT level and becomes small after 
GW correction. Above all, our scGWo predicts a wider 
range for M0S2 maintaining its direct gap characteristic 
compared to DFT results. This direct-indirect band gap 
transition also indicates that PL of monolayer M0S2 be- 
comes weaker under tensile tension of 1.41% and above, 
while corresponding experimental result is still lacking 
now. The optical gap shown in Fig. |1| will be discussed 
in the next subsection. 
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FIG. 4: Band gaps for monolayer M0S2 obtained by DFT, 
ScGWo, and BSE. Reported experimental (Expt.) 0, DFT 
[ p^[ , and scGW |^4j] results are also shown. The inset is the 
energy difference between the valence band tops at K and Y 
points. 



B. Excitonic effect in monolayer M0S2 

In this subsection, the optical properties of mono- 
layer M0S2 are discussed in details. From the techni- 
cal view, optical transition simulation needs the integra- 
tion over the irreducible Brillouin zone using sufficiently 
dense fc-point mesh. Naturally, the convergence of k- 
point sampling is important. First, for monolayer M0S2 
at strainless case (3.16 A), the optical adsorption spectra 
£2 (£xx =£yy) obtained by different fe-point meshes are il- 
lustrated in Fig. ||(a), in which the independent-particle 
(IP) picture is adopted within DFT (DFT-IP) and no lo- 
cal filed effect is included at hartree or DFT level. The 
first peak at about 1.78 eV is observed clearly in all the 
cases, corresponding to the K-K direct transition. The 
second significant peak located at about 2.75 eV is con- 
verged for 12 x 12 x 1 and 15x15x1 fc-point meshes. Other 
peaks in adsorption spectra between the two aforemen- 
tioned dominated peaks mainly originate from different 
irreducible k points with unequal weights in different k- 
point meshes. According to our analysis of projected 
density of states, the two significant peaks located in 1.78 
and 2.75 eV correspond to d-d and p-d transitions, re- 
spectively. Considering the dipolar selection rule only 
transitions with the difference Al = ±1 between the an- 
gular momentum quantum numbers I are allowed, i.e., 
the atomic d-d transition is forbidden. However, in the 
monolayer M0S2, due to the orbital hybridization, the 
VBM and conduction band minimum (CBM) still have 
p orbital contributions, especially the former; thus the 
VBM to CBM transition dominated by d-d transition is 
still allowed. As expected, the strength of this d-d tran- 
sition is weaker than the p-d transition as shown in Fig. 
|(a). 

As for the BSE calculations, in order to reduce the 
computational cost, we adopt 400 and 200 eV for the 
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FIG. 5: DFT-IP and scGW +BSE adsorption spectra for 
monolayer M0S2 at experimental lattice of 3.16 A (strainless 
case) obtained by different fc-point meshes. 



plane wave energy cutoff and response function energy 
cutoff (short for 400 and 200 eV for energy cutoffs), re- 
spectively, while the accuracy still can be guaranteed. 
Taking the strainless monolayer M0S2 for example, the 
scGWo band gap is 2.78 eV, resulting in only 0.02 eV 
difference compared to 2.80 eV aforementioned using 600 
and 300 eV for energy cutoffs. The calculated BSE spec- 
tra for strainless monolayer M0S2 are plotted in Fig. 
||(b). It is clearly that as fc-point mesh refines, the 
first peaks have a blueshift. For fc-point meshes 6x6x1, 
9x9x1, 12x12x1, and 15x15x1, the scGWq band gaps 
are 2.99, 2.84, 2.78, and 2.76 eV, respectively; the first 
adsorption peaks (optical band gaps) are 1.96, 2.08, 2.16, 
and 2.22 eV. Correspondingly, the exciton binding ener- 
gies are 1.03, 0.76, 0.62, and 0.54 eV, inferred from the 
difference between the QP (scGWb) and optical (scGWq- 
BSE) gaps. These calculated QP band gaps, optical 
gaps and exciton binding energies are also listed in Ta- 
ble ||. The convergence trend is obvious, particularly for 
the electronic band gap. However, due to the limitation 
of computation resource, scGWo calculations with more 
dense fc-point mesh are not performed here. Note that 
previous theoretical results showed a much larger value of 
exciton binding energy for monolayer M0S2. For exam- 
ple, a value of 0.9 eV for monolayer M0S2 (3.16 A) was 
obtained using empirical Mott-Wannier theory p3; and 
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TABLE I: QP band gap, optical band gap and exciton binding energy for monolayer M0S2 and WS2 are obtained by QP scGWo 
and BSE with and without spin-orbital coupling (SOC) adopting different energy cutoffs and fc-point mesh. All energies are in 
the unit of eV. 





Energy cutoffs 


fc-point 


E 9 


Eg (optical) 


Binding energy 


Monolayer M0S2 (3.160 A) 


400 and 200 


6x6xl(SOC) 


2.89 


1.87 


1.01 






6x6x1 


2.99 


1.96 


1.03 






9x9x1 


2.84 


2.08 


0.76 






12x12x1 


2.78 


2.16 


0.62 






15x15x1 


2.76 


2.22 


0.54 




600 and 300 


12x12x1 


2.80 


2.17 


0.63 


Monolayer M0S2 (3.190 A) 


600 and 300 


12x12x1 


2.66 


2.04 


0.62 


Monolayer WS 2 (3.155 A) 


400 and 200 


6x6xl(SOC) 


3.02 


1.97 


1.05 






6x6x1 


3.28 


2.21 


1.07 






9x9x1 


3.12 


2.34 


0.78 






12x12x1 


3.06 


2.43 


0.63 






15x15x1 


3.05 


2.51 


0.54 




600 and 300 


12x12x1 


3.11 


2.46 


0.65 


Monolayer WS 2 (3.190 A) 


600 and 300 


12x12x1 


2.92 


2.28 


0.64 



a value of 1.03 eV was obtained by Go Wb-BSE calcula- 
tions for monolayer M0S2 (3.18 A) using 6x6x1 fc-point 
mesh and including spin-orbital coupling p5f , which is 
the same as our above results using the same fc-point 
mesh without spin-orbital coupling. Apparently, these 
results overestimated the exciton binding energy due to 
the insufficient fc-point mesh in the first-principles calcu- 
lations. 

Experimentally, two close peaks observed in adsorp- 
tion spectrum of monolayer M0S2 around 1.9 eV are due 
to the valence band splitting caused by spin-orbital cou- 
pling. In our calculations, the spin-orbital coupling is 
omitted unless otherwise stated and this will not alter our 
main conclusions presented in the current study. In or- 
der to make a comparison, we also performed the scG Wb- 
BSE calculations with spin-orbital coupling using 6x6x1 
fc-point mesh and 400 and 200 eV for energy cutoffs. The 
two peaks in BSE adsorption spectrum locate at 1.87 and 
2.05 eV and the corresponding exciton binding energy is 
1.02 eV, consistent well with the aforementioned Go Wb- 
BSE calculations using the same fc-point mesh and en- 
ergy cutoffs while different pseudopotentials |2j| . Notice 
that the exciton binding energy obtained with and with- 
out spin-orbital coupling for monolayer M0S2 as shown 
in Table || is nearly the same, while the optical gap in the 
former case shifts about 0.1 eV towards lower energy due 
to the top valence band splitting of 0.17 eV according to 
our scGWo calculation. 

Based on our above results and discussions, the exci- 
ton binding energy for monolayer M0S2 is predicted to be 
in the range of 0.50—0.40 eV, while experimental value is 
still lacking. For one dimensional system like carbon nan- 
otube (CNT), a similar exciton binding energy of about 
0.4 eV was also obtained for semiconducting CNT with 



0.8 nm diameters (2^ ps| . 

For the evolution of exciton binding energy as a func- 
tion of strain, our results demonstrate that it is almost 
unchanged, i.e., 0.63 eV (strainless) , 0.62 eV (1% strain), 
0.62 eV (3% strain), and 0.59 eV (6% strain) (using 
600 and 300 eV for energy cutoffs and 12x12x1 fc-point 
mesh). The direct optical gaps are 2.17, 2.04, 1.81, 1.52 
eV for the four cases shown in Fig. g respectively. The 
former two are also shown in Fig. [| using the orange left 
triangles. The experimental optical gap for monolayer 
M0S2 was shown to be about 1.90 eV 0. Since there 
was no mention of specific lattice parameter, here it is as- 
sumed to be the strainless case as shown in Fig. |j. Notice 
that the consistency is good between our theoretical and 
experimental results. If spin orbital coupling is taken 
into account, the consistency will be improved further, 
since the first peak in the adsorption spectrum moves 
towards lower energy due to the top valence band split- 
ting. Most importantly, our results demonstrate that the 
optical gap of monolayer M0S2 is very sensitive to ten- 
sile strain, which can be tuned by depositing monolayer 
M0S2 on different substrates |L4j, whereas the exciton 
binding energy is insensitive to it. 



C. Chemical bonding properties of monolayer 
M0S2 

In order to gain further insight into the electronic 
structures under the tensile strain, we revisit the DOS 
shown in Fig. [j]. For the strainless case, the VBM states 
at K mainly originate from d xy and d x 2_ y 2, while under 
3% tensile strain, the VBM states at T mainly originate 
from d z 2 . The CBM at K is mainly contributed by d z 2 . 
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FIG. 6: MLWFs isosurface plots for Mo d z 2 (a and c) and 
S p z (b and d) for M0S2 under strainless (a and b) and 3% 
tensile strain (c and d) cases, respectively. 

The top valence bands at T point shift upwards compared 
to those at K point while CBM at K shifts downwards 
as the strain increases. Clearly under the tensile strain, 
the hybridization between Mo d and S p becomes weak. 
MLWFs can also illustrate the chemical bonding prop- 
erties of solids If]. MLWFs for Mo d z ? and S p z in 
monolayer M0S2 for both strainless and 3% tensile strain 
cases are plotted in Fig. ||. It is seen that the overlap 
of MLWFs for Mo d z i and S p z becomes weak due to 
the strain. In Fig. |(a), the MLWFs for Mo d z i have 
certain degree distortion towards S atom, while this is 
not observed in Fig. ||(c). As the strain increases, the 
ionic bond character is increased in monolayer MoS2- At 
about 10% tensile strain, according to DFT studies, a 
semiconductor to metal transition occurs |l3| . 

D. QP band structures and optical properties of 
strained monolayer WS2 

The QP band structures of monolayer WS2 under ten- 
sile strain are also investigated, motivated by its better 
performance than monolayer M0S2 used as a channel in 
transistor devices j3"0| . The calculation results are illus- 
trated in Fig. 0. Similar to monolayer M0S2, the scGWo 
QP band structures of monolayer WS2 also undergo a 
direct to indirect band gap transition as tensile strain in- 
creases. The direct band gaps for the strainless (at the 
experimental lattice of 3.155 A jl0|) and under 1% ten- 
sile strain cases are 3.11 and 2.92 cV, respectively, and 
the latter corresponds to the optimized lattice constant 
for monolayer WS2 from this work. The corresponding 
indirect band gaps under 3% and 6% tensile strains are 
2.49 and 1.78 eV, respectively. Note that for the strain- 
less case, our DFT result predicts monolayer WS2 to 
be an indirect band gap semiconductor with CBM only 
about 16 meV lower than the lowest conduction band 
at K points, which is contrary to recent full potential 
methods p0| . The difference may be originated from the 



(a)a=3.155A (b) a=3.190 A (c) a=3.250 A (d) a=3.344A 




FIG. 7: DFT, Go Wo, and scGWo QP band structures for WS2 
at lattice constants of (a) 3.155, (b) 3.190 (optimized lattice 
constant this work), (c) 3.250, and (d) 3.344 A , corresponding 
to 0%, 1%, 3%, and 6% tensile strain (with reference to 3.155 
A), respectively. The Fermi level is set to be zero. 

technical aspect of these calculations, such as the em- 
ployed pseudopotential method pi} ]. However, after the 
GW correction, a correct direct band gap is achieved. 

At the GW level, the dependence of energy difference 
between the VBT at K and T point on strain for mono- 
layer WS2 is similar to that of monolayer M0S2 as shown 
in Fig. [l| The direct to indirect transition for mono- 
layer WS2 is predicted to occur at 1.74% tensile strain 
by scGWq calculations. 

For optical properties of monolayer WS2, our calcu- 
lated QP band gaps, optical gaps and exciton binding 
energies are also listed in Table Q. It is obvious that the 
monolayer WS2 presents many similar properties com- 
pared to monolayer M0S2, for example, the gaps and 
exciton binding energy also demonstrate a convergence 
trend as fc-point mesh increases; the spin orbital cou- 
pling has little influence on the magnitude of the exciton 
binding energy. Notice that our scGWo calculation pre- 
dicts the top valence band splitting of monolayer WS2 to 
be 0.44 eV, larger than that of monolayer M0S2 of 0.17 
cV, because W is much heavier than Mo. The resulting 
first peak in BSE adsorption spectrum shifts 0.24 eV to- 
wards lower energy, also larger that that of monolayer 
M0S2 of 0.1 eV correspondingly. As for the strain effect, 
the BSE optical gap at our optimized lattice constant 
of 3.190 A is 2.28 eV, while at 3.16 A it is 2.46 eV, as 
shown in Table |[ The former corresponding to 1% tensile 
strain, results in 0.18 eV reduction of band gaps. This 
demonstrates that the band gaps and optical gaps are 
also very sensitive to tensile strain, whereas the exciton 
binding energy is not. Based on above analysis, we also 
predict the exciton binding energy for monolayer WS2 to 
be in the range of 0.5—0.4 eV, similar to that monolayer 
M0S2. Experimentally, the PL maximum of monolayer 
WS 2 locates between 1.94 and 1.99 eV [||. Considering 
the large shift of the peak in the BSE adsorption spec- 
trum caused by spin orbital coupling, our results at op- 
timized lattice of 3.190 A are consistent well with above 
experimental result. Notice that the self-consistent GW 
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TABLE II: Electron and hole effective masses (m*) derived from partially scGWo QP band structures for monolayer M0S2, WS2, 
MoSe2, and \VSe2 at different strains. The effective masses at K and V points are along KY and Mr directions, respectively. 







Compressive (1%). 


Experimental. 


Optimized. 


Tensile (3%) 


Tensile (6%) 


M0S2 


K, 


0.40 


0.36 (0.60 a ,0.35 6 ) 


0.32 


0.29 


0.27 




K h 


0.40 


0.39 (0.54 a ,0.44 6 ) 


0.37 








r, 








l.OO 


u.yu 


WS2 


K, 




0.27 


0.24 


0.22 


0.20 




K ; , 




0.32 


0.31 








r„ 








1.24 


0.79 


MoSe 2 


K e 




0.38 


0.36 








K h 




0.44 


0.42 






WSe 2 


K, 




0.29 


0.26 








K h 




0.34 


0.33 







"Effective masses listed here are averages of the longitudinal and transverse values in Ref. [E4j . 
Effective masses listed here are averages of the curvatures along the YK and KM directions in Ref. |2a 



method with vertex correction in W, which corresponds 
to the attractive interaction between hole and holes, also 
yields excellent band gaps pgfl . 

According to our above scGWo and BSE calculations 
for monolayer M0S2 and WS2, it is clear that the self 
energy within the scGWo calculations enlarges the band 
gap by accounting for the many body electron-electron 
interactions more accurately, while the strong excitonic 
effect results in a significant reduction of the band gap. 
Combining the two opposite effects on band gaps, the fi- 
nal resulting optical gap is consistent well with DFT band 
gaps. Therefore, the good band gap agreement between 
DFT and experiment is only a coincidence due to the fact 
that QP band gap correction is almost offset by exciton 
binding energy. This phenomenon was also observed in 
monolayer of hybridized garphenc and hexagonal boron 
nitride, which also have strong excitonic effect p4| . 

We also perform the scGWo QP band structures for 
monolayer M0S2 and WS2 under 1% compressive strains. 
Our results show that the compressed M0S2 has a direct 
band gap of 2.97 eV, while the compressed WS2 has an 
indirect band gap of 3.13 eV and K to K direct gap of 
3.30 eV. 

Our scGWo results show that both MoSe2 and WSC2 
are also a direct semiconductor at the strainless state. 
The experimental lattice constants Jl(J for MoSe2 and 
WSe2 are 3.299 and 3.286 A and the optimized lattice 
constants are 3.327 and 3.326 A, respectively; their di- 
rect K-K band gaps are 2.40 and 2.68 eV at experimen- 
tal lattices and 2.30 and 2.50 eV at the optimized lat- 
tices. Compared to the experimental lattice, the opti- 
mized lattice corresponds to 0.86% (1.22%) tensile strain 
for MoSe2 (WSe2), and the band gap also decreases with 
increasing tensile strain. 

E. Effective mass 

Based on the more accurate scGWq QP band struc- 
tures, the effective mass of carriers for TMDs are cal- 



culated by fitting the bands to a parabola according to 

fc2 j- 2 

E = 9 " n m , , where m e is the electron static mass. A 
fc-point spacing smaller than 0.03 A -1 is used to keep 
parabolic effects. Electron and hole effective masses (m*) 
at different strains arc collected in Table [fj. For MoS 2 un- 
der different strains, the CBM always locates in K point, 
and the electron effective mass K e increases with increas- 
ing compressive strain while decreases with increasing 
tensile strain. As for hole, initially the effective mass 
also decreases as the tensile strain increases. After the 
direct to indirect gap transition, VBM shifts to T with 
heavier hole, which also decreases as the tensile strain 
increases. Compared to the effective masses of 0.64 and 
0.48 for hole and electron at K point based on DFT cal- 
culation performed at the experimental lattice JTo[ for 
M0S2, the effective masses are reduced due to the GW 
correction in our study. For WS2, MoSe 2 , and WSC2, 
their masses also show similar behaviors. It is noted that 
at the same strain level, the mass of WS2 is the light- 
est, making WS2 more attractive for high performance 
electronic device applications since a lighter electron ef- 
fective mass can lead to a higher mobility. Theoretical 
device simulations also demonstrated that as a channel 
material, the performance of WS2 is superior to that of 
other TMDs Q. 



IV. SUMMARY 

In summary, the QP band structures of monolayer 
M0S2 and WS2 at both strainless and strained states 
have been studied systematically. The scGWo calcula- 
tions are found to be reliable for such calculations. Us- 
ing this approach, we find their band structures share 
many similar behaviors. At strainless state or under a 
small tensile strain, they maintain their direct band gap 
characteristics. As the tensile strain increases, they un- 
dergo a direct-indirect band gap transition. For the op- 
tical properties, exciton binding energy is predicted to 
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be in the range of 0.50—0.40 eV, which is much smaller 
than previous results and insensitive to the strains. Our 
calculated optical band gap is also consistent with ex- 
perimental results. In addition, we find that the electron 
effective masses of monolayer M0S2, WS2, MoSe2, and 
WSe2 reduce as the tensile strain increases, and WS2 
possesses the lightest mass among the four monolayer 
materials at the same strain. Importantly, the present 
work highlights a possible avenue to tune the electronic 
properties of monolayer TMDs using strain engineering 



for potential applications in high performance electronic 
devices. 
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